Information content of the low-energy electric dipole strength: correlation analysis 
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Background: Recent experiments on the electric dipole (El) polarizability in heavy nuclei have stimulated 
theoretical interest in the low-energy electric dipole strength, both isovector and isoscalar. 

Purpose: We study the information content carried by the electric dipole strength with respect to isovector 
and isoscalar indicators characterizing bulk nuclear matter and finite nuclei. To separate isoscalar and isovector 
modes, and low-energy strength and giant resonances, we analyze the El strength as a function of excitation 
energy E and momentum transfer q. 

Methods: We use the self-consistent nuclear density functional theory with Skyrme energy density functionals, 
augmented by the random phase approximation, to compute the El strength, and covariance analysis to assess 
correlations between observables. Calculations are performed for spherical, doubly-magic nuclei ^°^Pb and ^^^Sn. 

Results: We demonstrate that El transition densities in the low-energy region below the giant dipole resonance 
exhibit appreciable state dependence and multi-nodal structures, which are fingerprints of weak collectivity. The 
correlation between the accumulated low-energy strength and symmetry energy is weak, and dramatically depends 
on the energy cutoff assumed. On the other hand, a strong correlation is predicted between isovector indicators 
and the accumulated isovector strength at E around 20MeV and momentum transfer q ^ 0.65 fm~^. 

Conclusions: Momentum- and coordinate-space pattern of the low-energy dipole modes indicate a strong frag- 
mentation into individual particle-hole excitations. The global measure of low-energy dipole strength poorly 
correlates with the nuclear symmetry energy and other isovector characteristics. Consequently, our results do 
not support the suggestion that there exists a collective "pygmy dipole resonance," which is a strong indicator 
of nuclear isovector properties. By considering nonzero values of momentum transfer, one can isolate individual 
excitations and nicely separate low-energy excitations from the T — 1 and T = giant collective modes. That is, 
measurements at ^ > may serve as a tool to correlate the El strength with certain bulk observables, such as 
incompressibility or symmetry energy. 



PACS numbers: 21.10.Pc, 21.60.Jz, 21.65.Mn 24.30.Cz 

I. INTRODUCTION 

The electric dipole response of the atomic nucleus car- 
ries fundamental information on bulk nuclear proper- 
ties and shell structure [THS! and it also plays an im- 
portant role in nuclear reactions involving photo-nuclear 
processes. In particular, El strength impacts photo- 
absorption and radiative particle-capture processes oc- 
curring in stars during the cosmic nucleosynthesis [10] 
and the transmutation of radioactive nuclear waste pT] . 

The significance of the El strength, especially in the 
context of neutron-rich nuclei, has led to appreciable ex- 
perimental progress in this area p!2HT9] . Much excite- 
ment has been brought by a suggestion that the low- 
energy El strength, dubbed "Pygmy Dipole Resonance" 
(PDR) [20, "ST, is collective in nature, and can be under- 
stood as a motion of skin neutrons against the proton- 
neutron core. 

The current theoretical situation regarding the collec- 
tivity of the low-energy El strength is fairly confusing. 
While some papers advocate the existence of a collective 
PDR mode [T6l [TSl I22H25 j . having strength that is corre- 
lated with the nuclear symmetry energy fl6l |22l [23] , the 
collectivity of the low-energy El transitions, and their 
relevance to isovector nuclear matter properties (NMP), 



have been questioned by several other studies [26H3T] . 

In particular, the covariance analysis of Ref. [26^ has 
suggested the lack of correlation between PDR strength 
and nuclear isovector properties, such as neutron-skin 
thickness rgkin and symmetry energy. One reason for 
this was attributed to the local single-particle structure 
around the Fermi level (i.e., shell effects), which vary 
rapidly with global nuclear matter characteristics. On 
the other hand, the dipole polarizability - a global 
measure of the El strength - seems to be an excellent 
isovector indicator [I9l [26l [32l [33] . Similar conclusions 
have been reached in Ref. [28] , which concluded that the 
current data on the low-energy strength, and theoreti- 
cal predictions, are both too uncertain for a quantitative 
analysis. 

The aim of this work is to clarify the situation by 
means of the stringent correlation analysis based on 
methods of data analysis using least-squares optimiza- 
tion. We will consider, in particular, the correlations be- 
tween various nuclear observables (including NMP and 
properties of finite nuclei) and El strength. Our pa- 
per is organized as follows. Section [ll| describes our nu- 
clear EDF+RPA approach and motivates the choice of 
EDFs used. The definitions pertaining to El form fac- 
tors, strength functions, and resulting El sum rules are 
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discussed in Sec. |III[ The dependence of the low-energy 
El strength on selected nuclear matter properties is stud- 
ied in Sec. [Vl The results for transition densities are 



presented in Sec. IV, and Sec. VI contains the correla- 
tion analysis. While most calculations in this study were 
done for ^^^Pb, Sec. VII contains illustrative examples 
for ^^^Sn that demonstrate that our findings are general. 



Finally, Sec. VIII contains conclusions of the work. 



The focus of our survey is on properties of nuclear exci- 
tations in the electric dipole channel. To compute the El 
strength function, we employ the random phase approx- 
imation (RPA), using the same EDF as for the ground 
state to guarantee the full self-consistency. The details 
of RPA calculations follow Ref. 



III. STRENGTH FUNCTIONS 



II. MODEL 

The present investigation is based on the self-consistent 
nuclear density functional theory in the Skyrme-Hartree- 
Fock (SHE) variant [Ml [35]. The form of the SHE energy 
density functional (EDF) is derived from rather general 
arguments of a low- momentum expansion [36^ 37 . How- 
ever, since the coupling constants of EDF cannot be de- 
termined precisely from underlying nuclear forces, they 
are usually adjusted to selected nuclear data [34, 38, 39]. 
We use here two recent EDF parameterizations which 
try to embrace a large collection of various nuclear data 
[38l|40]. The survey [38] selected a pool of spherical nuclei 
which have been checked to be well described by a mean- 
field model [41 . The resulting parameterization SV-min 
was optimized by means of a least-squares procedure. Its 
(parabolic) least-squares landscape = X^{p)^ where p 
stands for the multitude of SHE parameters, carries im- 
portant information on the uncertainties, in particular 
on model features that are weakly determined by the op- 
timization process. This can be exploited by applying 
the standard covariance analysis as was done, e.g., in 
Refs. [26, 42 . To estimate the uncertainty related to the 
choice of fit observables, we alternatively look at the pa- 
rameterization UNEDEO [40 , which covers a much larger 
set of ground state data, also including many deformed 
nuclei. 

It is found that many aspects of EDF can be charac- 
terized by relatively few NMP, such as incompressibility 
isoscalar effective mass m* /m^ symmetry energy agym, 
and Thomas- Reiche-Kuhn (TRK) sum rule enhancement 
n (related to the isovector effective mass). By providing 
a set of parameterizations for which these NMP are sys- 
tematically varied [38j, one can track explicitly the in- 
fluence of the NMP on various nuclear properties. Base 
point of this set is the EDF called SV-bas, for which 
four NMP had been fixed: K = 234 MeV, m*/m = 0.9, 
<^sym = 30 MeV, and hz = 0.4 (chosen such that the giant 
dipole resonance (GDR) in ^^^Pb is described correctly). 
Starting from SV-bas, four sets of parameterizations were 
produced [38] by systematic variation of each one of the 
four NMP. It is to be noted that SV-bas and its system- 
atically varied sister parameterizations are adjusted with 
additional constraints on NMP. SV-min, on the other 
hand, is fitted only to data from finite nuclei and thus 
covers all uncertainty on NMP as left open by nuclear 
ground state data, a feature which make it ideally suited 
for our covariance analysis. 



In the following, we discuss various aspects of the elec- 
tric dipole strength, both isoscalar (T = 0) and isovector 
(T = 1). Detailed information on the strength is con- 
tained in the dipole transition form factor: 

A 

FP{q) = (^,1 Y,ji{qra)Y^o{n^)UT\^o) (1) 

where a runs over all nucleons, is the wave function of 
the z/-th RPA excitation, the ground state wave func- 
tion, q is momentum transfer, T is the isospin quantum 
number, and IIt stands for the isospin projector: 

N ^ Z - 

no = Ilprot+nneut 7 ni = ^Ilprot ~ ^Ilneut- (2) 

(For spherical nuclei, such as ^^^Pb and ^^^Sn considered 
in this work, it is sufficient to consider the single mode 
with the magnetic quantum number ja = 0.) The q- 
dependence carries the pattern of an excitation mode. 
A complementing coordinate-space view is provided by 
radial transition density, which is the Fourier transform 

POO 

) (r ) = 4^ / dr j, {qr) (q) . (3) 

Since the RPA excitation spectrum of a heavy nucleus 
such as ^^^Pb is fairly dense, it is more convenient to 
consider the energy distribution of form factor strength: 

SPiq,E) = ^Y1 \FP{Q)?Gr{E - E,), (4) 

where E^, is the RPA excitation energy of 
state u and Qr{F — E^) is a Gaussian fold- 
ing function having energy-dependent width 
r{E^) = max[0.2MeV,(^^-8MeV)/6MeV] - to 
simulate nuclear damping effects, which increase with 
energy. We note small F (and improved resolution) in 
the low energy region, which is the very focus of this 
work. 

The factor q~'^ in Eq. Q is introduced to cancel 

the leading g -dependence of Fj) \q) at small momen- 
tum transfer. Indeed, at g ^ 0, the isovector strength 
S^p\q^E) is proportional to that of the familiar dipole 
operator D = rYioIIi, which we denote by Sd{E). The 
next-order term in q leads to the "compressional dipole" 
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operator D^^"^ = (r^ — ^{r'^)r)YioflT, the leading ob- 
servable in the isoscalar channel [45H49] . and the corre- 
sponding strength is denoted by Sjj (E). For energy cuts 
at finite we shall consider the general dipole operator 
F = ji{qr)YiotlT as in Eq. (B. 



UNEDFO 



dipole T=1 

compress. T=0 

connpress.T=1 



c 

c 
q 

T3 



C 





15 20 25 30 35 
E (MeV) 

FIG. 1. (Color online) Isoscalar and isovector El strength 



Sd (solid line), (dotted line), and (dashed line) in 
2°^Pb predicted with (a) UNEDFO, (b) SV-bas, and (c) SV- 
min EDFs. The region of GDR is marked by vertical dot- 
ted lines; these limits are somewhat arbitrary as the GDR 
strength does not show clear energy bounds. 

To assess the importance of the low-energy El 
strength, we introduce the accumulated strength: 



A^{E)= dE' E'" SoiE'), 
Jo 



(5) 



where O can mean either the mere dipole compres- 
sional dipole Dc^ or general dipole F at a given value of 
q. In Eq. ([5|, we use a constant folding width F = 0.5 
MeV in So- For the dipole operator, the inverse-energy- 
weighted sum rule A^i{oo) is related to the electric dipole 
polarizability [5 



(6) 



ueRPA 



while the energy- weighted sum rule AP{oo) is the TRK 
sum rule. 

Figure [l] shows El strength Sd{E) and S^^^ (E) in 
208pb predicted with UNEDFO, SV-bas, and 'sV-min 



EDFs. The isovector dipole strength Sd is concentrated 
in the GDR region of 12-16 MeV. The systematic shift 
of the GDR strength between EDFs refiects the change 
in TRK sum rule enhancement (/>:=0.4 for SV-bas, 0.25 
for UNEDFO, and 0.08 for SV-min) [38 . The isovector 

compressional strength 6*^^ roughly follows Sd- It does 
also contain a high energy branch above 30 MeV which, 
however, is not of interest in the context of this study. 
The isoscalar compressional strength is fairly broad, with 
a low-energy component at 5-18 MeV and a high-energy 
concentration at 20-30 MeV [8 . 

Using the standard representation of El strength as 
shown in Fig. [T] it is difficult to separate the low- 
energy strength from the tails of giant resonance modes. 
A better separation can be obtained by considering 
nonzero values of i.e., by analyzing the full distri- 
bution sf\q,E). Fi sure 2 shows the structure of the 
El strength Q - together with individual proton and 
neutron strength distributions - for SV-bas in the {E^ q) 
plane. As expected, the isoscalar channel shown in 
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FIG. 2. (Color onhne) El strength Q for ^°^Pb computed 
with SV-bas as a function of E and q. (a) Neutron strength; 
(b) proton strength; (c) T = 1 strength (scaled by 2); and (d) 
T — strength (scaled by 0.25). Since Sf rapidly decreases 
with ^, rather than the actual strength Q we show the q- 
scaled strength S^P {q, E) = {q, E){1 + aqq^f with aq = 
9fm^ 

Fig. [2|d) has no strength at g' ^ 0, as the spurious 
isoscalar collective center-of-mass El mode is removed 
in our calculations. The large isoscalar strength ap- 
pears at somewhat larger momentum transfer around 
q = 0.65fm^. Therein, one can clearly see the two 
branches of the isoscalar dipole compressional mode at 6- 
13 MeV and 24 MeV. The GDR centered around 13 MeV 
dominates the T = 1 strength at low values of It seems 
that the group of excitations between 6 and 11 MeV has 
predominantly isoscalar compressional character. They 
seem to be higher-g extensions of some low-energy isovec- 
tor dipole strength (often referred to as FDR strength). 
These low-energy dipole transitions can thus be consid- 
ered as shadows of the low-energy dipole compressional 
modes. The neutron and proton El strengths displayed 
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in Figs. [2|a) and (b), respectively, show that the states 
having large transition form factors are all mixed proton 
and neutron excitations. While the neutrons carry more 
strength generally, pure neutron modes do not appear. 

The upper branch of the compressional mode around 
22-26 MeV shown in Fig. |2] is not so interesting for our 
study. Therefore, in the following, we will restrict the 
energy range. Figure [3] illustrates the low-energy pattern 
of the El strength for three EDFs of Fig. [l] It is en- 
couraging to see that the general behavior of El strength 
is EDF-independent, so the detailed discussion of RPA 
modes will be based on SV-bas calculations. In the en- 
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FIG. 3. (Color online) Similar as in Fig. |2] for T = (left) 
and T = 1 (right) El strength computed with (a) UNEDFO, 
(b) SV-min, and (c) SV-bas EDFs. 



tive to study this point in the coordinate-space in terms 
of the transition density pi^'' of Eq. In order to 

present results as a function of we performed a Gaus- 
sian folding as in Eq. Q. Figure [i] depicts the general 
behavior of El transition density in^^Pb as a function of 
E and r. At £; < 12 MeV, p^^^ varies rapidly with E, ex- 
hibiting a complex multi-nodal behavior in both isospin 
channels (cf. discussion in Refs. [33l [50]). The strong 
state dependence, together with the lack of a common 
pattern, are both fingerprints of weak collectivity. A dif- 
ferent situation is seen in the GDR region. Although the 
energy band of 12-17 MeV covers a large number of RPA 
states, the T = 1 transition density has the same pattern 
throughout, with a pronounced surface bump. A simi- 
lar pattern is seen for the T = compressional mode at 
22-28 MeV. Such coherent behavior is characteristic of a 
collective mode. 
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ergy range below 10 MeV, only a handful of modes can be 
seen. Actually, we find more RPA states (e.g., 20 modes 
between 6 and 10 MeV) in that range, but only a third 
of them are visible in El excitation. These strong-El, 
low-£^ modes have a fairly diverse character. The low- 
est states below 7.5 MeV are predominantly isoscalar in 
character, having some contribution to the T = 1 chan- 
nel. The mode at 8 MeV has a mixed character (isovector 
at low g'- values and isoscalar at higher q), while the mode 
at 9 MeV is predominantly an isovector excitation. The 
states around 10 MeV are again dominated by T = 0. 
Their T = 1 component is weak, and it is mainly concen- 
trated at larger i.e., these states do not carry significant 
dipole {q = 0) strength. 

IV. TRANSITION DENSITIES 

The interesting question is to what extent the lowest 
El modes can be considered as collective, as discussed 
(or assumed) in some studies. To this end, it is instruc- 



FIG. 4. (Color online) Energy- averaged El transition density 
in ^°^Pb calculated with SV-bas as a function of E and r for 
T = (top) and T = 1 (bottom). 



V. LOW-ENERGY ELECTRIC DIPOLE 
STRENGTH 

Although transition density and form factor reveal a 
world of local details, it is desirable to quantify their 
information content in terms of a few integrated quan- 
tities. One such quantity is the accumulated low-energy 
El strength, 

B{El;E,n..)= B{El,v), (7) 

which is proportional to the accumulated strength 
^o*(^max)- This observable is somewhat ambiguous as 
it depends on an assumed cutoff energy £^max- Figure [5] 
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shows the predicted B{E1; E^g,^) in ^^^Pb as a function 
of two isoscalar (incompressibihty K and isoscalar effec- 
tive mass m^/m) and two isovector (symmetry energy 
ttgym and TRK sum rule enhancement k,) NMP indicators 
[26^ . The calculations were carried out with four sets of 
EDF parameterizations around SV-bas [38] by systemat- 
ically varying the NMP of interest. We used the value 
^max=10MeV that corresponds to the beginning of the 
GDR region in SV-bas, see Fig. [3) 
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FIG. 5. (Color online) Dependence of the low-energy El 
strength ^ on selected NMP {K,m* /m, asym, i^) in ^°^Pb 
predicted with SV-bas family of EDFs. The cutoff energy is 
^max=10MeV. 



When considering each variation independently, one 
can see a monotonic trend in each case. But does it mean 
that a well-defined correlation exists between the low- 
energy El strength and the bulk properties considered? 
Clearly, the answer to this question cannot be given based 
on Fig. [5] alone, as B{E1] £^max) seems to depend on sev- 
eral NMP, both isoscalar and isovector in character, and 
the plot says nothing about the possible coupling be- 
tween them as an underlying EDF parametrization is 
systematically varied. Second, if a parametrization is too 
constrained in a given interaction channel, all variations 
probing this channel are correlated by construction. (An 
example of such a situation was given in Ref. f26^ in the 
context of a relativistic mean field model RMF-(5-t, too 
constrained in the isovector channel to be used in a mean- 
ingful correlation analysis.) Although the explicit track- 
ing with respect to single NMP, as done in Fig. [5j can 
be very useful to uncover hidden dependencies between 
observables, it does not explore all conceivable variations 
and is thus insufficient to quantify true correlations. A 
more exhaustive method will be discussed in the next 
section. 



VI. CORRELATIONS 

Thus far we have analyzed the structure of dipole 
modes in terms of transition form factor and transition 
density. These structures do not show convincing signa- 
tures of collectivity in the low-energy region below GDR. 
As a more global measure, we inspected the integrated 
low-energy dipole strength, which indeed displays some 
complex dependence on NMP. For a more exhaustive 
(and quantitative) analysis, we will now exploit the strin- 
gent method of covariance analysis and scrutinize the cor- 
relations between integrated dipole observables and col- 
lective bulk properties (NMP, dipole polarizability, and 
neutron skin). 

As already mentioned in section the empirical ad- 
justment of EDF parameters produces a nearly-parabolic 
landscape near the minimum, i.e., x^(p) ^ Xo + 
"^Zij PiPj^PiPjX^ ' The vicinity of the minimum in the 
p-space, for which x^(p) = Xo + considered as 

the space of acceptable parameter variations (la region). 
Observables are also functions of the EDF parameters: 
A = A{p). A correlation between observables A and B 
within a given model can then be assessed by means of 
the correlation coefficient 

Cab \ (8) 
V AA2 AB^ 

where the overline means an average over the space of ac- 
ceptable parameters [24l [26l [42] . Such correlation analy- 
sis had been proven useful in previous surveys. For exam- 
ple, it was concluded that the electric dipole polarizabil- 
ity is a good isovector indicator that strongly correlates 
with the neutron radius of ^^^Pb [24, 26j. 

The correlation between El strengths in ^^^Pb and ba- 
sic nuclear matter properties (K, m*/m, agym, is dis- 
cussed in Figs. [6]|8] as a function of E. Figure |6] il- 
lustrates the correlation with the accumulated dipole 
strength A^{E) The case of n = —1 is shown in 

Fig.mFa). As noted earlier, A^i{E) is the accumulator for 
the dipole polarizability a^- It is interesting to see how 
the sensitivity to agym (and insensitivity to other NMP) 
develops at high energies. At GDR energies and below, 
however, the correlations vary dramatically. The changes 
are particularly rapid around the lower end of the GDR 
region (9-10 MeV), which is often used as a cutoff energy 
^max to determine the low-energy El strength ([7|. It 
is apparent that even small changes in £^max may change 
correlations. This result casts serious doubts on standard 
(cutoff-dependent) definitions of the pygmy strength. 

The corollary of the emerging specific sensitivity is seen 
in Fig.|6|b), which illustrates the case of n = 1. At large 
energies, the accumulated dipole strength Ai{E) devel- 
ops into an unambiguous measure of k,. As expected, K 
and m* /m never correlate with A-^^{E)^ regardless of the 
energy region. This is to be expected as these NMP are 
isoscalar indicators [26 . 

The summed dipole strength Aq (E) is often used as a 
measure of the pygmy mode. We show in Fig. [7] how 
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FIG. 6. (Color online) Correlation between the energy- 
weighted accumulated strength (E) ^ in ^°^Pb and se- 
lected NMP (incompressibility isoscalar effective mass 
rrf /m^ symmetry energy cisym, and TRK sum rule enhance- 
ment Hi) as a function of E obtained with SV-min. Top: 
n — —1; bottom: n — 1. The GDR region is marked by 
vertical dotted lines. 



it correlates with NMP. At zero momentum transfer, 
Fig. [7|^a) 5 the correlation pattern strongly resembles that 
of Fig. [6|b) , including strong correlation with k at high 
energies. The three cases displayed in Figs. |6]and[7|^a) 
have certain common features. Namely, they all show (i) 
rapid changes near the lower and the upper end of the 
GDR region; (ii) strong correlation with n in the GDR 
region; and (iii) medium correlation with ttgym 

around 

E = 7 MeV. 

A stronger handle on agym is provided by the summed 
dipole strength Aq (E) at nonzero momentum transfer. 
Figure illustrates the case of integrated form factor 
strength dit q = 0.65 fm~^. A particularly large correla- 
tion with agym is predicted in a wide energy range above 
the GDR region and below 25 MeV. A strong correlation 
with hi in the GDR region still holds at this value of q. 

The complementing n = isoscalar strengths are 
shown in Fig. [sj^a). At low values of the summed 
compressional dipole strength Aq''{E) contains little in- 
formation on NMP. The sensitivity to agym at the low- 
E region, suggested in [33 , is minor. Figure [sj^b) 
demonstrates that some sensitivity to K appears at 
q = 0.65 fm~^, especially at low energies below the GDR 
region and in the high energy region above 20 MeV. 

We now turn to a correlation between isovector form 
factor strength Aq{E) and the four strong isovector in- 
dicators: the slope of symmetry energy ftgy^^ at the satu- 
ration density, slope of the neutron EoS Eneut/N at neu- 
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FIG. 7. (Color online) Similar to Fig. [g] but for the accu- 
mulated dipole strength Aq(E) for two values of momentum 
transfer: (a) q = and (b) q = 0.65 fm~^. 
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FIG. 8. (Color online) Similar as in Fig.|7|but for the T = 
compressional dipole El strength. 



tron density 0.08 fm^. El polarizability q^d, and neutron- 
skin thickness rgkin in ^^^Pb. Recall that these four ob- 
servables are highly correlated with each other and with 
agym [26 . Consequently, all four correlations displayed in 
Fig. |9] are practically identical, and also agree with the 
results for agym shown in Fig. [7| the summed isovector 
El strength above the GDR at the intermediate values of 
momentum transfer around q = 0.65 fm~^ is an excellent 
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isovector indicator. As illustrated in Fig. |9]jc), however, 
this correlation deteriorates at higher values of q. 
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FIG. 9. (Color online) Correlation between the T — 1 accu- 
mulated El strength in ^"^^Pb and selected isovector indica- 
tors (ttgym, ^neut/A^, Q^D, Tskin) in SV-miu as a function of E 
for ^ = (a); q = 0.65 fm"^ (b); and q = LMfm"^ (c). 



Fig. 10 illustrates the information content of the ac- 
cumulated inverse-energy-weig hted El strength A^-^{E) 
with respect to the four isovector observables discussed in 
the context of Fig. [9] Again, we see that the correlations 
are very similar for all four observables and the T = 1 
accumulated strength at g = 0.65 fm~^ nicely correlates 
with isovector indicators around E = 2{) MeV. The accu- 
mulated inverse-energy-weighted isoscalar strength does 
not seem to correlate well with isovector observables. 



VII. A QUICK GLANCE AT ^^^Sn 

It is interesting to check whether the findings obtained 
for ^^^Pb are specific to this nucleus or whether they are 
of a more general nature. To this end, we inspect the case 
of a neutron-rich, doubly-magic nucleus ^^^Sn. Figure 11 
shows its T = 1 and T = El strength predicted with 
SV-bas. The overall pattern, with GDR transitions in 
the isovector channel and the dominant dipole compres- 
sional mode in the isoscalar channel, closely resembles 
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FIG. 10. (Color online) Correlation between the T — 1 
(left) and T = (right) accumulated inverse-energy-weighted 
El strength A^i in ^"^^Pb and selected isovector indicators 
(ttgym, ^neut/A^, ttD, Tskin) in SV-mlu as a function of E for 
q — (top) and q — 0.65 fm""*^ (bottom). 
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FIG. 11. (Color online) Similar as in Fig. [2] panels (c) and 
(d), except for ^^^Sn. 



that shown in Fig [2] for ^^^Pb. In particular, the low- 
energy modes are isoscalar at \ow-q values and isoscalar 
at higher momentum transfer. 

Figure 12 shows the transition densities for ^^^Sn. 
Again, the pattern is extremely similar to the case of 
^^^Pb. Finally Fig. 13 shows the correlations of the four 
basic NMP with accumulated isovector dipole strength. 
It is, again, similar to the analogous Fig. l7[a) for ^^^Pb. 
We also carried out calculations for other neutron-rich 
nuclei, such as ^^^Sn, and our conclusion is that the case 
of ^^^Pb is representative to other spherical closed-shell 
heavy systems. 



VIII. CONCLUSIONS 

To clarify the ambiguous theoretical situation with 
regard to the interpretation of the low-energy electric 
dipole strength, we study its information content with 
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FIG. 12. (Color online) Similar as in Fig.ji] except for ^^^Sn. 
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FIG. 13. (Color online) Similar as in Fig. |7|a), except for 
^^^Sn. 



respect to isovector and isoscalar indicators characteriz- 
ing bulk nuclear matter and finite nuclei. As a theoretical 
tool, we use self-consistent nuclear density functional the- 
ory augmented by the random phase approximation. The 
inter-observable correlations are computed by means of 
covariance analysis. We use well-calibrated Skyrme en- 
ergy density functionals that allow for systematic varia- 
tions of isoscalar and isovector parameters characterizing 
nuclear matter properties. (The extension of the correla- 
tion analysis to other models, including relativistic EDF, 
will be carried out in a forthcoming study [51, •) To avoid 
ambiguities due to strong shell effects and pairing, our in- 
vestigations have been limited to the doubly-magic heavy 
nuclei ^^^Pb and on ^^^Sn, and the results were very close 
in both cases. To separate various electric dipole modes, 
we study the El strength as a function of excitation en- 
ergy E and momentum transfer q. By going to the q- 



dimension, we are able to resolve excitation modes that 
are close in energy, and can nicely separate low-energy 
excitations from the T = 1 and T = giant resonances. 

Our study fully confirms the main conclusion of our 
previous work [26j. Namely, the low-energy dipole ex- 
citations cannot be interpreted in terms of a collective 
pygmy dipole resonance mode associated with the mo- 
tion of skin neutrons. The detailed inspection of transi- 
tion form factor and density shows a pattern of rapidly 
varying particle-hole excitations. Moreover, it seems that 
a large fraction of the low lying states are isovector "shad- 
ows" of the underlying low-energy part of the isoscalar 
compressional mode. 

At the level of more global observables, we have studied 
the El strength integrated up to a given cutoff energy. As 
a first indicator of correlations, we have plotted in Fig. [5] 
the changes of the low-energy strength with systemati- 
cally varied nuclear matter properties. This has led to 
a complex picture where every NMP has some influence 
on the low energy strength. To assess correlations in a 
meaningful way, we went beyond the simple inspection 
of trends and employed the standard tools of statistical 
analysis. While the local El strength at low energies 
weakly correlates with isovector indicators such as sym- 
metry energy, the magnitude of those correlations varies 
rapidly around the lower end of the GDR region. This 
strong dependence on the cutoff energy indicates that a 
pygmy strength cannot be unambiguously defined; hence, 
it is not a particularly useful quantity. 

Of particular importance are the fully integrated 
strength; the much celebrated sum rules. The T = 1 sum 
rules, such as the dipole polarizability (the inverse 
energy weighted sum rule), the sum rule enhancement 
factor n (strongly correlated with the energy weighted 
sum rule), and the total accumulated El strength are 
all excellent probes of nuclear isovector properties. At 
the isoscalar side, we found a relatively strong corre- 
lation between nuclear incompressibility and the accu- 
mulated isoscalar dipole strength for q = 0.65 fm~^ (at 
E = 9 MeV and above 20 MeV). In general, experimental 
studies at > in the isovector and isoscalar channel 
could provide an excellent window for a better charac- 
terization of El strength and relating it to fundamental 
nuclear properties. 
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